Diluted magnetic semiconductor quantum dots: an extreme sensitivity of the hole 
Zeeman splitting on the aspect ratio of the confining potential 
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The valence band states confined in infinitely deep quantum dots made of diluted magnetic semi- 
conductors (DMS) are considered theoretically. A complex anisotropic structure of the valence bands 
in DMSs with cubic symmetry described by the full Luttinger Hamiltonian is taken into account. It 
is found that the Zeeman splitting is very sensitive to the shape of the confining potential and, in 
particular, to its orientation relative to the direction of an external magnetic field. This sensitivity 
has its origin in a mixing of different spin components of a hole wave function which takes place for 
finite hole wave vectors k. Several consequences of the effect are discussed, including a possibility to 
control the inter-dot tunneling by an external magnetic field. It is shown also that the polarizations 
of optical transitions in a single DMS quantum dot depend on details of geometry of its confining 
potential as well as on the strength of the magnetic field. 



I. INTRODUCTION 

From the moment of initial stages of their investiga- 
tion in mid-seventies (see, e.g., Jl[, |[ [|) until present 
diluted magnetic semiconductors (DMSs) and quantum 
structures involving them, continue to attract a consid- 
erable attention of the scientific community (see, e.g., 
H U). The interest has recently gained impetus with 
new focus on spintronic materials (see, e.g., recent re- 
view ||). The part of this interest stems from the fact 
that the presence of an exchange interaction between the 
band carriers and electrons from the partially-filled d- or 
/-shells of magnetic ions in DMSs gives rise to a unique 
possibility to control electrical and optical properties of 
quantum structures by means of an external magnetic 
field in a range difficult to achieve in other materials. 

This s-d exchange interaction, in the mean field ap- 
proximation, can be cast pj into Zeeman-like form 



iJ s _ d = N Q ax(S) • s, 



(1) 



in the case of the conduction electrons, and the p-d in- 
teraction in the same approximation 
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p—d 



N Px(S) ■ J, 



(2) 



in the case of the valence band holes. Here § and J are 
operators of the electron spin and the total angular mo- 
mentum of the hole, respectively, x is the molar frac- 
tion of magnetic ions introduced substitutionally into the 
semiconductor matrix and (S) represents their average 
localized spin. 

The corresponding exchange constants Nqo. for the 
conduction band states and N f3 for the valence band 
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states in majority of II- VI DMSs are usually of op- 
posite signs @. Usually also, the absolute value of 
Nq(3 is considerably greater than Nga. For example, in 
Cd!_ x Mn x Te N a = 0.22 eV and N P = -0.88 eV [§. 
Thus the giant spin splitting (GSS) of excitonic states in 
DMSs is mainly determined by the spin splitting of the 
hole states. 

The excitonic GSS in DMS quantum structures is de- 
termined by the hole splitting to an even greater extent 
than in bulk crystals. The reason for this is that the 
electron and the hole states in low-dimensional semicon- 
ductor structures are energy dependent mixtures (stem- 
ming from the fc-dependence of the mixing in the bulk) 
of the Bloch amplitudes with Tq and Fg symmetries of 
the center of the Brillouin zone point group representa- 
tions. Since the negative exchange constant Nq/3 is sev- 
eral times greater in magnitude than the positive N a, 
even a small admixture of Tg states to Tq state can re- 
duce significantly (up to 30% in the case of quantum 
wells) the value of total effective exchange constant of 
the conduction electrons |J. For quantum wires (QWRs) 
and, especially, for quantum dots (QDs) this reduction of 
the conduction band exchange constant is expected to be 
even more significant. A degree of a reciprocal process 
of a reduction of the value of the effective exchange con- 
stant of the holes due to an admixture of Tq symmetry 
to Tg states is less pronounced because of the mentioned 
difference of the magnitudes of the exchange constants 
Noa and Nq(3. Therefore, the exciton spin splitting in 
DMS quantum structures, which combines the splitting 
of the electron and the hole comprising it, is dominated 
by what happens in the valence band. 

An example of such processes, that affects the Zeeman 
splitting of the holes, is an additional heavy hole-light 
hole mixing induced by the QD confining potential. This 
has been noted and exploited in Ref. BO] where author 
used an exactly solvable model, namely of the Luttinger 
Hamiltonian in a spherical approximation with spatial 
confinement given by a spherically symmetric quantum 
dot potential. The reduction of the Zeeman splitting 
of an optically active ground state of an exciton due to 
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heavy-light hole mixing was found to be independent of 
the QD size. The degree of the reduction depends only 
on the relation between the heavy and the light hole ef- 
fective masses in the crystal and for CdTc it amounts to 
about p w 0.8. 

In the present paper we pursue the study of the con- 
sequences of the heavy hole-light hole mixing induced by 
the localization of the holes by a QD potential. However, 
we concentrate mainly on the dependence of the mixing 
on the shape or, more specifically, on the aspect ratio of 
the confining potential. We demonstrate that in different 
QDs, giving rise to similar confinement energies, the sym- 
metry and, thus, the total angular momentum of the hole 
ground state might be completely different. To illustrate 
this point, let us imagine that a characteristic length of 
a quantum dot in, say, z-direction L z is much smaller 
than that in x- and y-directions L z -C L x , L y . The situ- 
ation becomes then similar to the case of a quantum well 
(QW) , where the ground state of the hole is purely of a 
heavy hole in character with J z = ±3/2. In the oppo- 
site case, when L z 3> L x ,L y the structure resembles a 
QWR with the ground state of holes being mainly the 
light hole in character due to the greater masses in x-y 
plane jflL |l2). Yet, from the point of view of the confine- 
ment energy only both these two extreme cases may be 
the same. 

One of the consequences of different symmetries of the 
hole ground states in different structures having the same 
confinement energy is a strong dependence of the expec- 
tation value (J z ) on the shape of the confining poten- 
tial. On the other hand, the Zeeman splitting of the hole 
states produced by the interaction Hamiltonian (|2|) in the 
first order of perturbation theory is proportional to the 
value of ( J z ) (we shall always assume that external mag- 
netic field B that induces the magnetization appearing 
in Eqs. (0) and (d) is along the z-axis). One can ex- 
pect, thus, in the case of DMSs a strong dependence of 
the effective exciton spin splitting on the aspect ratio of 
quantum dots. We will show also that due to the strong 
anisotropy of the hole states the very orientation of the 
QD with respect to an external magnetic field affects dra- 
matically the Zeeman splitting. 

It was shown in Ref. |l3| that due to the ^-dependence 
of s(p)-d exchange constants their values in quantum 
structures are reduced to a certain extent. Neverthe- 
less, for the sake of simplicity , to describe s-d and p-d 
exchange interactions we use in this paper the values of 
these constants equal to their values in the bulk. This 
simplification becomes invalid only in very small QDs 
(with radius smaller than about 30A), while here we fo- 
cus our attention on much greater structures and we show 
that even in those the dependence of the Zeeman split- 
ting on their geometry is significant. Admittedly, for 
small QDs one should replace bulk No (3 by the values of 
p-d exchange constants, calculated in [ 
change, of course, the qualitative resu 
present paper. 

To sum up, we are considering here the shape-induced 
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changes of ( j z ), and their various consequences, in DMS 
QDs while neglecting the variation of exchange constant 
Nq(3 itself since we expect it to be of smaller importance. 



II. THEORETICAL MODEL 

We consider the valence band states in a diluted mag- 
netic semiconductor QD in a shape of a rectangular par- 
allelepiped with dimensions L x , L y , and L z . This partic- 
ular choice of the form of confining potential is motivated 
by two reasons. First, by changing the relations between 
L x , L y and L z one can control the asymmetry of the con- 
fining potential. This is convenient both for calculations 
and for interpretation of the results obtained. The second 
reason is that the symmetry of QD potential in this par- 
ticular case is rather low. In a general case, this is D^h 
symmetry. Such a low symmetry prevents appearance 
of additional constants of motion, such as the angular 
momentum which "creeps in" in the case of spherically 
symmetric QDs. 

For simplicity we consider the case of a dot with in- 
finitely high potential barriers surrounding it. The pa- 
rameters L x , L y and L z can be treated then as char- 
acteristic length scales of the hole localization in a real 
QD (i.e., with finite potential barriers) in x, y and in- 
directions. 

Our basic starting Hamiltonian is the Luttinger Hamil- 
tonian jL4| describing the hole states in a cubic semicon- 
ductors which can be presented in the form 
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273 K X Ky(J X Jy 



JyJx 



(3) 

k x k z (j x j z + j z j x ) 



K z Ky (J z Jy ~t" JyJz) 



where K x = —id/dx, K y = —id/dy, K z = —id/dz, with 
the z-axis is parallel to [001] crystallographic direction, 7$ 
are the Luttinger parameters, while J x , J y and J z being 
4x4 operator matrices of the projections of the angular 
momentum J = 3/2 in the basis of Tg Bloch amplitudes 
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(5) 
(6) 
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We expand the orbital part of each of the spin com- 
ponents (Q)-(j^) of the hole wave function in the basis of 
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the solutions of the Schrodinger equation for an infinitely 
deep QD. The hole wave function is then represented by 



* = 



+3/2 N x N y N. 

^ EEEd 



L x L y L z 



J z = -3/2n=lm=l 1=1 

x sin.(nk w x) sin(mk y y) sixi(lk z z) ■ ijij m , (8) 

where k x = ir/L x , k y = n/L y , k z = n/L z , N x , N y and 
N z are (finite) numbers of the basis functions used in the 
expansion and A"' m ' 1 being numerical coefficients. Let 
us note that in general is not an eigenfunction of J z or 
any other component of J. 

For the case of the hole ground state let us first limit 
the expansion (||) to the terms with m = n = I — 1. 
In this case the orbital part of the hole envelope has the 
form 



L x L y L- 



■ sin(fc x x) sin(fe^y) sin(fc z z). (9) 



This approximation is equivalent to neglecting of term 
proportional to 73 in the Hamiltonian (||) or, more 
strictly, to a reduction of the spherical term oc 73 (K • J) 2 
to a one with cubic symmetry oc 73 (K 2 J 2 + K 2 J 2 + 

K 2 J 2 ). In fact, the terms proportional to K 2 in Eq.(||) 
do not mix states with different orbital functions in the 
expansion (||) and the hole ground state has the form (0). 

With the above approximation the last term in Eq. (B|) , 
being linear in the operators K a , has only vanishing 
matrix elements. At the same time, the second term 
in Eq.(^) is proportional to the squares of the matri- 
ces J a . Thus, the non- vanishing matrix elements of this 
term involve states with AJ Z — ±2 only. As a result, 
the 4x4 matrix of the Hamiltonian (0) splits into two 
2x2 identical submatrices in the basis (+3/2, —1/2) and 
(—3/2, +1/2). These matrices have the form 



H = 



P + Q R 
2m V R P- 



(10) 



where 



Q = 72 (k 2 x + k 2 y 
R = -V3 l2 (kl - k 2 y ) 



kl) 
2k 2 ) 



In the case of a strictly cubic QD (k x =k y =k z ), the 
ground state of the hole is four-fold degenerate. The 
hole states can be characterized by the projection of the 
angular momentum J z = ±3/2, ±1/2. We can remove 
this degeneracy by reducing the symmetry of the Hamil- 
tonian from Oh to Dih (k z ^ k x = k y ). This leads to the 
heavy hole-light hole splitting. J z remains a good quan- 
tum number. The hole ground state doublet is then the 
heavy hole doublet with J z = ±3/2 if k^+k 2 < 2k 2 , while 
it corresponds to the light hole doublet with J z — +1/2, 
if k 2 x + kl > 2k 2 . 



In the general case, k x ^ k y ^ k z , (the symmetry of the 
Hamiltonian being Dih) the hole states are the mixtures 
of the components with J z — ±3/2 and J z = +1/2. In 
the absence of an external magnetic field the ground state 
with the energy 



E = 



2m 



P- 



VQ 2 + R 2 ) 



(ii) 



is twofold degenerate. This effective spin doublet is char- 
acterized by an expectation value of the projection of the 
total angular momentum 



(J, 



± 



,1 + (d+v^rry 



where 
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V3|fc 2 



(12) 



(13) 



For degenerate states the choice of the wave functions is 
not unique. To obtain Eq.(|l2"|) we use the correct zero- 
order functions with respect to the perturbation GJ Z . 

For each value of k x and k y it is possible to find a value 
of k z that leads to vanishing of (J z ). The corresponding 
condition is 5 = l/y/3. In the case of k x > k y to satisfy 
this equality one has to choose k z — k y . This result easy 
to see: the equality k z = k y under the condition k x > k y 
means that the hole ground state is a superposition of 
the states with J x — ±3/2 only. The average value (J z ) 
for such states is equal to zero, similarly to the case when 
the states are the mixtures of states with J z — ±3/2 and 
then (J x ) and (J y ) vanish. 

It is seen, thus, that the expectation value of J z for 
the hole ground state in a QD varies in a wide range 
from (J z ) — ±3/2 to (J z ) — depending on the exact 
shape of the confining potential, i.e., relation between its 
extension in three crystallographic directions. 

The situation changes somewhat if we take into ac- 
count more than one term in the expansion (||). The 
projection J z ceases to be a good quantum number for 
the holes even in the case when k x = k y . This is due 
to the presence of the term proportional to 73 in Eq. (||) . 
Moreover, the hole states now are the mixtures of all four 
spin components. Nevertheless, as before, one can vary 
effectively the value of ( J z ) in the ground hole state from 
nearly ±3/2 to by changing the shape of the confin- 
ing potential. To illustrate this feature we plot in Fig. [j] 
the expectation value of the projection of the total an- 
gular momentum J z for the ground spin doublet in the 
QD having L x — 50A and L y = 70A (solid line) and 
L x = 40A and L y = 80A (dashed line) as a function of 
L z . Hereafter we use the number of basis functions in the 
expansion (JsJ) N x = N y = N z =6, which were checked to 
be sufficient to ensure numerical accuracy (about 1 meV 
in calculation of hole energy). In our calculations we use 
the CdTe Luttinger parameters 71 = 5.3, 72 = 1.62 and 
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30 40 50 60 70 80 90 

L [A] 

FIG. 1: The expectation value of the projection of the total 
angular momentum (J z ) as a function of the QD size in the 
z-direction. The solid line corresponds to the structure with 
L x — 50A and L y = 70A, while the dashed one corresponds 
to the structure with L x — 40A and L y = 80A. 

73 = 2.1 p5[ . Figure [l] shows that the value (J z ) is, 
in fact, very sensitive to the form of the quantum dot 
potential. 
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FIG. 2: The Zeeman splitting of the hole ground states 
in Cdo.97Mno.03Te DMS quantum dots in the presence of a 
magnetic field B || z at T = 2K. The solid lines correspond 
to the structure with L x = 50A, L y — 70A and L z — 65A, 
the dashed lines - to the structure with L x — 65A, L y — 70A 
and L z = 50A. 

As mentioned above, the value of the Zeeman splitting 
of the valence band states in the magnetic field B || z is 
proportional, in the first order of perturbation, to (J 2 ). 
However, the perturbation theory is not quite justified 
if the value of the Zeeman splitting becomes compara- 
ble to the energy separation between the size-quantized 
states in a QD. To calculate the Zeeman splitting in such 
case we perform a numerical diagonalization of the full 
Hamiltonian matrix for the holes, which is represented 



by the sum of the Luttinger Hamiltonian and the p-d 
exchange interaction (Eh. Note, that we neglect here the 
direct influence of the magnetic field on the hole states. 
The magnetic field in our model just induces the mag- 
netization on the sample, which is proportional, in turn, 
to the average spin (S) of the magnetic ions appearing 
in (Q). In the calculations that are presented below we 
use the parameterization of the value (S) in real DMS 
structures, given in [p~G|] . 

III. DISCUSSION 

The evolution of the positions of the two lowest energy 
levels in the valence band in Cdo.97Mno.03Te DMS quan- 
tum dots having L x = 5(M, L y = 70A and L z = 65A 
(solid lines) and L x = 65A, L y = 70A and L z = 50l 
(dashed lines) in the magnetic field B || z is presented in 
Fig. |^. Such QDs provide exciton confinement energy of 
about 400 meV, that is of the correct magnitude as com- 
pared to existing experimental observations in real self- 
assembled QDs made of this material [jiTj. In a small 
magnetic field the value of the Zeeman splitting of the 
ground state is proportional to the expectation value of 
the projection of the total angular momentum of the hole 
on the direction of the magnetic field. This quantity is 
(J z ) = 0.003 for the first structure and (J z ) = 1.2 for 
the second structure, which is just the first one merely 
rotated by 7r/2 around the y-axis. Due to the fact that J z 
is not a good quantum number, the magnetic field B || z 
induces the mixing of the hole states corresponding to 
different size-quantized levels. With an increase of the 
magnetic field the admixture of the higher energy states 
to the ground levels becomes significant and it produces 
a deviation of the value of the Zeeman splitting from that 
predicted by the first order of perturbation theory. Nev- 
ertheless, Fig. H shows that due to a strong dependence of 
(J z ) on the actual geometry it is possible to obtain DMS 
quantum dots with the same, homogeneous distribution 
of the magnetic impurities, that have the same space- 
quantized energies but, nevertheless, differ essentially in 
their effective g-factors. 

This sensitivity of the effective g-factor to the exact 
shape and orientation of the confining potential can be 
used to control hole tunneling between DMS quantum 
dots. To clarify this point, let us consider two QDs lo- 
calizing spatially the holes in DMS QW in two in-plane 
directions. In the third direction, say ^-direction, the 
holes are localized due to the presence of a QW confin- 
ing potential. Such a pair of the QDs is, thus, made 
of the same material and have the same dimensions in 
the x-direction. The remaining two spatial dimensions 
of these QDs are chosen in such a way as to result in 
different expectation values of some component of the 
angular momentum of a hole, say (J z ), in the plane of 
the QW. Moreover, the ground state energy of the hole 
in the QD having the greater value (J z ) should be some- 
what higher than that in the other QD, see Fig. ^. The 
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FIG. 3: Schematic view and the evolution of the ground 
energy of the hole states in two Cdo.97Mno.03Te quantum dots 
in a magnetic field B || z. Label 1 corresponds to 70A x 
6uA x 5uA quantum dot, label 2 to 70A x 50A x 65A quantum 
dot. The resonant field, that couples the two quantum dots 
enabling the hole resonant tunneling, is about 1.2 Tesla at 
T=2K. 



energy difference between these two ground states and 
the distance between the QDs may prevent the tunnel- 
ing of a hole between these two QDs. Therefore, we refer 
to these QDs as uncoupled in this case. Now, the mag- 
netic field B || z shifts the ground states of holes in QDs. 
This shift is different for each dot and is proportional to 
( J z ). This is due to the fact that the ( J z ) are different for 
our two dots. At some magnetic field energies of ground 
states of holes in two QDs may become the same, see 
Fig. §. One can achieve, thus, the condition favorable 
for the resonant tunneling. Let us assume now that the 
hole originally is localized in the first QD (e.g., where the 
ground state in the absence of the field is lower) and that 
the time length of the magnetic field pulse is made equal 
to one half (three halves, five halves etc.) of the period of 
the hole oscillations between the two QDs. Then, after 
such a pulse the hole will be localized in the second QD. 

The system is scalable to some extent. It is possible, 
for example, to add a third QD with another value and/or 
direction of the resonant magnetic field that couples the 
second and the third QD by the tunneling. As a result, 
by applying the magnetic field that effectively couples the 
second and the third QD one can also transfer the hole 
between these dots without affecting the hole present in 
the first QD. Such an arrangement can be envisaged as 
useful in quantum logic gates. 

Of course, the quantitative analysis of a possibility of 
such hole tunnel transport as well as the quantitative de- 
scription of the Zeeman splitting of excitons in DMS QDs 
requires the consideration of the finite potential barriers 
and the valence band-conduction band mixing. 

Along with the Zeeman splitting, the polarization 
characteristics of optical transitions in QDs have to be 
strongly affected by the shape of the confining potential. 
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FIG. 4: The probability of optical transitions associated with 
the fourfold degenerate exciton ground state as a function of 
L z for the QD with L x = 7uA and L y = 50A. The three lines 
correspond to three different linear polarizations of the light. 

One can expect this even by inspecting Fig. |l|. A strong 
dependence of the character of the hole ground state (i.e., 
the relative contribution of different spin components to 
the wave function) on the shape of a QD should result 
not only in a variation of (J z ) but also of probabilities 
of optical transitions. The values describing the prob- 
abilities of various optical transitions in different linear 
polarizations of light associated with the fourfold degen- 
erate exciton ground state as a function of L z for QD 
having L x = 70A and L y = 5f)A is plotted in Fig. |J. 
With an increase of L z , the situation changes from QW- 
like (L z < L x ,L y ) to QWR-like (L z > L x ,L y ). At the 
same time, the intensity of the optical transition in the 
E || z polarization increases, while the intensities of the 
remaining two polarizations, considered here, decrease. 

Finally, let us mention yet one more effect that takes 
place in DMS quantum dots. It was shown in Ref.|j"^| 
that the localization of the holes in two spatial directions 
in the case of DMS quantum wires results in a possibility 
to control the character of the ground state in the valence 
band by means of a magnetic field. This leads, in turn, to 
a strong magnetic field dependence of the polarization of 
optical transitions in such structures. A similar situation 
occurs, of course, in the case of DMS quantum dots. In 
Fig. H we plot the magnetic field dependence of the prob- 
abilities of optical transitions in circular er + , cr~, and 
in two linear polarizations for ground (non-degenerate) 
exciton state in Cdo.97Mno.03Te QD with L x = 70A, 
L y = 60A and L z = 50A. It is seen, that by appli- 
cation of a magnetic field one can effectively control the 
polarization characteristics of the light emitted from such 
DMS QD. 

Let us emphasize, that we use in our calculations the 
full, anisotropic Luttinger Hamiltonian and a rectangular 
profile of a QD, so the selection rules for the optical tran- 
sitions differ substantially from those obtained assuming 
spherically symmetric quantum dots. For example, in 
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FIG. 5: The probability of optical transitions in two circular, 
and two linear polarizations for the ground state of exciton 
in Cdo.97Mno.03Te DMS QD with L x = 70,4, L y = 60A and 
L z = 50A in the magnetic field B || z at T=2K. 

our case the ground state of the exciton in a QD is not 
a dark state. Moreover, all valence band states in our 
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fects are much more pronounced, as a direct consequence 
of the large values of the effective g-factor of the holes. 
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